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FIRST PASSAGE TIMES FOR COMBINATIONS OF RANDOM LOADS 



by 

P.A. Jacobs 

Operations Research Department 
Naval Postgraduate School 
Monterey, CA 93943 

0 . ABSTRACT 

Structures are subject to changing loads from various sources. 

In many instances these loads fluctuate in time in an apparently 
random fashion. Models are considered for which the stress put 
on the structure by various loads simultaneously can be described 
by a regenerative process. The distribution of the first time 
until the stress on the structure exceeds a given level x is studied. 
Asymptotic properties of the distribution are given for a large 
stress level x and for the tail of the distribution for fixed 
finite stress level x. Simulation results are given to assess 
the accuracy of using the asymptotic results to approximate the 



distribution . 



1 . 



INTRODUCTION 



Many physical structures are subject to varying physical 
loads from various sources: wind, snow, and earthquakes being 

examples. In many instances the total load experienced by a 
structure varies in time in an apparently random fashion. Certain 
load components, e.g. those of snow and ice accumulation, vary 
rather slowly; others such as those associated with winds or 
earthquakes occur more nearly as impulses. The problem is to 
design structures to withstand a coincidence of such loads with 
approximately a prescribed (high) probability. In engineering 
terms we wish to work towards developing a rational safety factor 
criterion for designing structures to withstand the combination 
of loads anticipated. 

In this paper we will study the time for the load combination 
process, {X(t); t >_ 0 } to exceed a given stress level for the case 
in which the load combination process is a regenerative process 
taking non-negative values (c.f. Cinlar [1975], page 298). Many 
of the models for load combinations that have been studied are 
regenerative processes; (cf. Pearce and Wen [1983], Wen [1981], 
Shanthikumar and Sumita [1983]) . 

EXAMPLE 1. The load combination process is the superposition 
of two load types: shock loads and constant loads. A constant 

load (e.g. caused by snow or rain accumulation) exists at one 
level for a random time and then changes to a new level. Let 
Y (t) denote the constant load magnitude at time t. The instants 
of change in the magnitude of Y(t) occur according to 
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a renewal process with inter-renewal distribution H; successive 
magnitudes of the constant load process are independent iden- 
tically distributed with distribution F. Impulse or "shock" 
loads (e.g. caused by wind gusts or earthquakes) occur at 
random moments and have time duration of length 0. Let Z(t) 
denote the shock load magnitude at time t. Given the constant 
load magnitude at time t, Y(t) = y, the probability a shock 
load will occur in the time interval [t,t+h] is y(y)h + o(h); 
the magnitude of the shock load is conditionally independent 
of everything given Y(t) = y and has distribution G(y;*). 

The load combination process magnitude at time t is 
X(t) = Y ( t ) + Z (t) . { X ( t ) ; t >_ 0} is a regenerative process 

with regeneration times the times of change of the constant 
load process. 

Gaver and Jacobs [1981] studied a special case of the 
above model in which the interrenewal distribution H is exponen- 
tial, and the shock load process is a compound Poisson process 
independent of the constant load process with shock arrival 
rate y and magnitude distribution function G. Other similar 
load combination models have been studied in the past; cf. 

Peir and Cornell [1973], Wen [1977], Pearce and Wen [1983]. 

Let X ( t ) denote the magnitude of the load combination 
process at time t. The process {X(t) , t _> 0 } is assumed to 
be a regenerative process. Let 

T = inf { t > 0: X(t) > x} , (1.1) 

the first time the load combination exceeds a level x. 
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Section 2 is concerned with asymptotic distribution of T 
as x -* oo. It is shown that under certain assumptions the distri- 
bution of the normalized random variable T (E[T ] is approximately 
unit exponential for large x, and error bounds on the rate of 
convergence are obtained. This result is related to that of 
Keilson [1979, page 134]. 

In Section 3 the tail of the distribution of T for finite x, 

x 

p{T >t}, will be studied. An asymptotic result concerning the 
exponentiality of P{T x > t} for large t is given. 

In Section 4 simulation results are presented to study the 
accuracy of using the two . asymptotic results to approximate the 
distribution of T x . 

2. THE DISTRIBUTION OF THE FIRST PASSAGE TIME 
FOR THE LOAD COMBINATION PROCESS 

Let { X ( t ) ; t ^0} be a regenerative process taking non-negative 

values representing the load combination process; [cf. Cinlar [1975], 

page 298] . Let S^ denote the nth regeneration time; (S n ) is a 

renewal process. We will assume there is a regeneration at time 

0. Let T be as in (1.1) and put <j>U) = E[exp(i£T )], the Fourier 
X x 

transform of T . A renewal theoretic argument yields 

<MO = E[exp{i£T x );T x < Sj + E [ exp { i£T x } ; S 1 < T x ] (2.1) 

= E [exp{ i£T } ; T < S,] + E [exp{ i£S, } ; S, < THU) 

XX-*- _L _L X 



Thus , 
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<M5> 



( 2 . 2 ) 



i£T 

E [e X ;T x < S 1 ] 

ns: 

1 - E[e ■ L ;S 1 < T x ] 



Similarly, if m(x) = E[T ], then 



m (x) = E [T ; T < S, ] + E[T ; T > S, ] 
x x — 1 x x 1 



(2.3) 



E [T ; T < S, ] + E[S, ;T > S, ] 
x x — 1 l x 1 



+ P{ T > S. }m(x) 
X -L 



Therefore , 



E[min(T ,S )] E [min (T , S. ) ] 

m(x) = 1 - PIT > S, } = PIT < S. } ' (2.4) 

x 1 x — 1 

The following assumptions will be made for the remainder 
of this paper. 



0 < E[S^] < oo ; (2.5) 

lim P { T < S. } = 0 . (2.6) 

x — 1 



It now follows from (2.4)- (2.6) that 

lim m (x) = °° . (2.7) 

x-«o 



4 



LEMMA 1. 



TTf 1 * ( srnrr - rvcd E[s i :S i < V 

E[sh 

± e(x) 2E [min( Sl ,T x )] + E[S 1 ] 

3/2 E[S l} 1/2 

+ 7 wr -- -to- T 



lJ • —Whu ) 1 ' 2 



2E [min (S. ,T ) ] 

-L X 



where B(x) = P{T < S, } • 

Proof. It follows from (2.2) that there are random variables 
and such that 



<f>( 



m(x) 



E [1 + 



iC 

m(x) 



T + 

x 



c • T < 

e l'x - 



( 2 . 8 ) 



[1 - E [ 1 + 



iC 

m ( x) 



+ £ , 



Vi 



-l 



= [P{T x 



S- ) 



+ 



i 5 

m (x) 



E [T 



x 



x 



< S.^ 



+ E[e 1 ;T x < S x ] ] 



x [P{T < S. } 

X _L 



iC 

m(x) 



ElS^S-j^ 



< T ] 
x 



E [e 2 ;S 1 < T x ] ] 



= (m(x)P(T < S, } + i£E(T ; T £ S, ) + m(x)E[e,;T £ S ]) 

X X X -L X X 

X [m(x)P{T x £ S x } -iCE[S 1 ;S 1 < T x ] - m(x)E[e 2 ;S 1 > T^] ] 
= (a + i£b) x (c - i£d) 1 . 
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The difference between the Fourier transform of T ,(E[T ] 1 ) 

x x 

and that for a unit exponential is 
[(a+Ub) x (c-iCd)" 1 ]- (1-iC)" 1 (2.9) 

= [(a-c+^ 2 b) + i£(b-a+d)] x [ (c-i£d) ( l-i£) ] 1 

_ NUM 
" DEN ‘ 

The term 

a -c +C 2 b = m(x)P{T x £ S ± } + m(x)E[e 1 ;T £ S ± ] (2.10) 

- m(x)P{T < S , } + m ( x ) E ( e 0 ; S , < T ] 

x — 1 2 1 x 

+ C 2 E [T ; T < Sj 
x x — 1 

= m(x)([E[e 1 ;T x £ S ± ] + EU^-S-j. < T^] ) 

+ ^ 2 E[T ;T < SJ. 

XX -L 

It follows from (2.4) that the term 

b-a+d = E[T ; T £S,]-[m(x)P{T <S,}+m(x)E[e 1 ;T < S, ] ] 

X X -L X J. IX 1 

+ E [ S 1 ; S 1 < T x ] (2.11) 
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b-a+d = E[min(T ,S )] - E[min(T ,S,)] - m(x)E[e.,;T < Sj 

X X _L x — _L 



= -m(x)E[e^;T x £ . 



Thus , 



NUM| = (m(x) (E[e 1 ;T x £ + E [ e 2 ; S ± < T] ) 



(2.12) 



+ £ 2 E[T x ;T x £ s ] _] } 2 + C 2 (-m(x)E[e i; T x < S 1 ] ) 2 • 



From Lemma 1 on page 512 of Feller [1971] 



Et I e i i ; T x i s i J i l ( ffrur ) 2 E[T x :T x i s i ] 



(2.13) 



E[ I £ 2 i ;S 1 ' T x' i 5<5 (It )2e[s i ;S i < V 



(2.14) 



Thus 



l» UM l 5 i<sfer> (EtT x :T x i s i ] + Els i ;s 1 < T x' > 



(2.15) 



+ C E[T x ;T x < S 1 ] + U | 



mfxT I E|T x :T x i Si 1 * 



Further , 



c - i£d| > |C| E[S 1 ;S 1 - T x ] 



(2.16) 



7 



Thus 



a + i£b 
c 



- r^I?l E[s r s i < V 



(2.17) 



i sUf I lE[T x ;T x i S l> + E|S 1' S 1 < T x>> 



+ 1*1 E[T x ;T x EIT X ;T X i s l>- 



Application of Schwartz' Inequality yields 



ITT ‘ r^ir 



E[S, ; S , < T ] 



l'~l 



x 



(2.18) 



< E[S?] + (1 + olfL ) E [T„;T„ < Sj 



— 2m(x) 1- 



2m (x) x' x — 1 



1 -"' 2 ’ ■ ■ _ U _ L ) ( e [ S 2 ] 1//2 P { T < S, } 1//2 ) . 



— 2m(x) 1 



E[Sf] + (1 + 



2m(x) 



'x - 1 



Let 6(x) = 



P{T £ S } . Then 
X 1 



m(x) = 



E[min(T x ,S 1 ) ] 

(Tlx) 



(2.19) 



Thus , 



TTT 1 ^ ( srlo ) ■ T^rr 1 E[S 1 ;S 1 < T x ] 



( 2 . 20 ) 



< B(x) 



+ id 



B(x) 3/2 e[s 2 ] 1/2 



2E [min (S . , T ) ] 

-i- X 
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THEOREM 1. Under Assumptions (2.5)— (2.6) 



l P{ HTnb" - y} ' (1 -e Y ) I < 0(P{T X £ S 1 } 1/1 °) . 



Proof . It follows from Lemma 1 that 

T 

Ul" 1 e lC | E[exp{iC ^^ ^| ^ }] - (1-U)" 1 ) |E[S 1 ;S 1 < T x ] 

1 iet ^i^vv] + E ^ 1/2 e(*> 1/2 ) 

2 3/2 E ts ^] 1/2 

+ IM 2 (6(x) 3 / — ) 

X X 

where B(x) = P(T £ S n } as before. 

X X 

Thus, applying Lemma 2 on page 538 of Feller [1971] for 
T = B(x) -1/5 , 



tt|p{T _< y } “ (l-e y )| <_24B(x)'*" // ^ 

X 



+ E[S 1 ;S 1 < T x HB(x) 2/5 [B(x) 



E[S 1 ] 2 1/2 
2E [min ( S . ,T ) ] + E[S 1 ] 

X X 



e[s 2 ] 1/2 

+ b ( x) -3/5 1 e (x) 3/2 — X ] 

+ B (xj 3 BIX) 2E [min (S . ,T ) ] J 

X X 



The result now follows. 



. 21 ) 



B (x) 1/2 ] 
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EXAMPLE. In Example 1 



x oo _ 

8(x) = F ( x) + / F(dy) / H (dt ) [ 1 - exp{-y (y ) G (y , x-y ) t } ] 

0 0 



x 

= 1 - / F(dy)h(y(y)G(y, x-y) ) (2.22) 

0 

- 00 -st 

where h(s) = j e H(dt) is the Laplace transform of H and 
0 

F ( t ) = 1 - F ( t ) . If C and S are two independent random variables 
having distributions F and H respectively, then 

3(x) = 1 - E [exp {- y ( C) G ( C , x-C) S }; C £ x] . (2.23) 

Consider two independent load combination processes 
{ ( t ) ; t _> 0} and {X 2 (t);t _> 0} of the type in Example 1. 

Assume that the conditional distribution of the shock load process, 
given the constant load process, is the same for both load combinati 
processes; the constant load magnitudes have distributions F^ 
and F 2 respectively with F-^(t) £ F 2 (t); and the times between 
constant load changes have distributions and respectively 
with H^(t) <_ ^(t) . Let and have distributions F^ and 

1L respectively for i = 1,2. If y ^ y(y)G(y,x-y) is an 
increasing function of y for y _< x, then 

-y (C )G(C , x-C. )S 

8 1 (x) = 1 - E [e . ;C 1 < x] 
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6 1 (x) < 1 



-y (C-)G(C ,x-C 0 )S 0 
- E[e 2 2 2 ;C 2 < x] 

5 B 2 (x) . (2.24) 

Thus, Theorem 1 suggests that the convergence to exponential 
of the distribution of the time for the load combination process 
to exceed a level x is faster for load combination process 1 
than that for process 2. This behavior has been seen in the 
simulation studies reported in Section 4. 

3. THE TAIL OF THE DISTRIBUTION OF T FOR FINITE x 

x 

In this section we will study the behavior of P{T > t} 

2 x 

for finite x and large t. The probability, P{T^ > t} , 
satisfies the following renewal equation 

t 

P{T > t} = P{S n > t,T > t} + f L (du) P{T > t-u} (3.1) 
x l x r , J x x 



where 



L (t) = P{S, < t , T > S, } . (3.2) 

x 1 — x 1 

We will assume L x (0) = 0. 

Following the argument on page 376 of Feller [1971] we 
will assume that for each x there exists a constant k (x) such 
that 
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00 

/ e K(x)u L (du ) - 1 , (3.3) 

0 x 



Y (x) = / te K{x)t L (dt) (3.4) 

0 X 



is finite, and the function 



g (t) = e K(x)t P(S > t , T > t} (3.5) 

X _L X 



is directly Riemann integpable. 

It now follows from the Key Renewal Theorem that 



lim e K{x)t P{T > t} = / g ( t ) dt . (3.6) 

t-M» X ' 0 X 



Since g is assumed to be directly Riemann integrable, 
integrating (3.3) by parts yields 



oo 

1 = P { T > S. } + k ( x) / e K(x)t P{ S. > t , T > S }dt • (3.7) 

X _L q -L X X 



Thus 



P{T < S ) 

k(x) = — ^ 1 . (3.8) 

I e k ( x ) t p ^ s > t>T > s } dt 

n 1X1 
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It follows from (2.4) that 



k (x) m (x) 



E [min ( T , S. ) ] 



(3.9) 



oo 



0 



/ e K (x) 1 P{S. > t,T > S. }dt 
; lxl 



The defining equation for k(x), (3.3), (3.8), and (2.6) 

imply that as a function of x, k(x) is nonnegative and decreasing 

with lim k(x) =0. If it is further assumed that there 

x->°° 0S^ 

exists 0 > 0 such that E[e 1 ] < °°, then it follows from the 
dominated convergence theorem and (3.9) that 



EXAMPLE. In Example 1, assume that H(t) = e , 
p(y) E p, and G(y,x) E G(x). Then k (x) satisfies the equation 



lim k (x)m(x) = 1 . 

X-M50 



(3.10) 



-At 



x 



A 



1=J F(dy) 
0 



A + pG (x-y ) - K (x) 



(3.11) 



where G(x) = 1 - G(x). 



Further 



, r X < (x) t - t - G(x-y)t . 

g^(t) = J e e e F(dy) . 



0 



Thus 
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00 X -j 

/ g (t)dt = / F(dy) — 

0 0 A + yG(x-y) - k (x) 



Therefore , 



lim e K< ' x)t P { T > t} 

. x 

-£->00 



_1 

Ay (x) 



where 



°o — 

y ( x) = / / F (dy ) t e K(x) Ae" At e _yG(x_y)t dt 

0 0 



X 



= / F(dy) 



A 



(A + yG (x-y) - k (x)) 



2 * 



If it is further assumed that A = y = 1, then it 
from (3.9) and (3.11) that 



k ( x) m (x) = 



/ F(dy) -± 

0 1 + G ( x-y ) 

X 

/ F (dy) = 

0 1 + G ( x-y ) -k (x) 



X 



/ F (dy ) = 

0 1 + G (x-y) 



Thus k(x) < — — r in this case. 
- m(x) 



(3.12) 



(3.13) 



(3.14) 



follows 



(3.15) 
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In the further special case in which F(x) = G(x) = e x , A = 1 
and y = 1, k(x) satisfies the equation 



x 



1 = f(K(x)) 



>-y 



1 + e 



- (x-y ) 



- k(x) 



dy 



(3.15) 



where 

f ( 1 ) = |[e X - e" X ] (3.16) 

— x 

and for y ^ 1 such that 0 . <_ y £ 1 - e , 

f ( y ) ( 1-y ) 2 = (1-y) ( l-e“ X ) -e _X [ In ( ( 1-y ) +e -X ) -In ( ( 1-y ) e _X +e~ X ) ] . 

(3.17) 



Further , 



y (x: 



= / 

0 



x 



(1 + e 



-(x-y) 



(x))‘ 



dy ; 



or 



Y (x) (1-k (x) ) 3 = ( 1-k ( x) ) ( 1-e X )-2e x [ In ( ( 1 -k (x) ) +e X ) 
- In ( (1-k (x) ) e x +e x 



)] • 



(3.18) 



4. A SIMULATION STUDY OF THE ACCURACY OF THE EXPONENTIAL 
APPROXIMATIONS TO THE FIRST PASSAGE TIME DISTRIBUTION 

In this section some results of a simulation study of the 

accuracy of using the asymptotic results of Sections 2 and 3 

to approximate the distribution of T are reported. All 

simulations were carried out on an IBM 3033 computer at the 

Naval Postgraduate School using the LLRANDOM II random number 

generating package (see Lewis and Uribe (1981)). 

The model of Example 1 was simulated for various cases 

of distributions H, F, G and shock arrival rate y . Each 

realization simulated the sample path of the process (X(t); 

t > 0} and the first passage times 0<_T £ T <_..._< T 

12 n 

were recorded for several levels x. < x~ < . . . < x . The number 

1 — 2 — — n 

of replications was 5,000. Sample moments and quantiles were 

computed for each T . A more detailed account of the simula- 

i 

tion can be found in Noh (19 84) . 

4.1 The Exponential Approximation P{T > t} = exp{-E[T ]t} 

X X 

a. Model A 

Tables 1-1, 1-2 and 1-3 report the simulation results 
for the model in which H = F = G are all unit exponential and 
y = 1 . Table 1-1 reports simulated sample mean and coefficient 
of variation for T^ for various levels of x. As expected, 
the simulated coefficient of variation approaches the exponen- 
tial distribution's value of 1 as x gets larger. To assess 
the quality of the exponential approximation 1 - exp{E[T^] "'"t} 
to the distribution of T , quantiles from the simulated data 
were computed and compared to the approximating exponential 
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TABLE 1-1 



Simulated Moments for T in the Case F = G = H = exp(l) and y = 1 



x-Level 


E[T ] Coeff Var [T ] 


0.50 


0.276 1.984 

*(0.01) 


1.00 


0.669 1.451 

(0.014) 


2.00 


2.000 1.177 

(0.033) 


3.00 


4.878 ' 1.080 

(0.074) 


4.00 


11.539 1.024 

(0.167) 


5.00 


26.945 1.002 

(0.382) 



* 

( ) 


is the standard error of the mean 



17 



quantile 



q = -E [ T ] ln(l-a) 

a x 



(4.1) 



where for this model E[T ] can be computed analytically and is 

X 

found in Gaver and Jacobs [1981]. Table 1-2 reports the 
simulated quantiles with the approximating exponential quantiles 
below in parentheses. As expected, the exponential approxi- 
mation is better for the large level x = 5 than for x = 0.5. 

The approximation is also better for a > 0.5. 

One way the distribution of differs from an exponential 
is that it has an atom at 0 ; in particular for the model A 
p (t = 0} = 1 - F(x) = e x . The sample quantiles for the 

X 

simulated conditional distribution of T given T >0 appear 
in Table 1-3. Below the simulated quantiles in parentheses 
appear approximating exponential ones computed as 



The quantiles of the simulated conditional distribution 
are much closer to their exponential approximation than those 
for the unconditional distribution, 
b. Model B 



the model which is the same as Model A except that the distri- 
bution of the constant load magnitudes F is exponential with 
mean 1/2. The simulated means and coefficient of variations 






a 



c 




(4.2) 



Tables 2-1, 2-2, and 2-3 report simulation results for 
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Simulated a-Quantiles for the Conditional Distribution of T x , Given 
T > 0, for the Model with F = G = H = Unit Exponential and y = 1 
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the approximating exponential quantile 



TABLE 2-1 



Level x 
0.5 

1.0 

2.0 

3.0 

4.0 

5.0 



Simulated Moments 



E [T ] 
x 



0.577 
* ( 0 . 010 ) 



1.329 

( 0 . 022 ) 



3.980 

(0.058) 



10 .530 
(0.150) 



28.014 

(0.395) 



75.425 

(1.070) 



for Model B 

A 

Coeff Var [T ] 
x 

1.453 
1.144 
1 .036 
1.006 
0.997 
1.003 



*( ) is the standard error of the sample mean 



Simulated a-Quantiles for the Distribution of T for Model 
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is the approximating exponential quantile 
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Approximating exponential quantiles 



appear in Table 2-1. An analytical expression for E[T ] for 
this model appears in Gaver and Jacobs [1981]. This expression 
is used in the approximating quantiles (4.1) and (4.2). The 
approximating quantiles for the conditional distribution of T 

X 

given T >0 which appear in Table 2-3 are closer to their 
corresponding simulated quantiles than the approximating quan- 
tiles for the unconditional distribution. Comparison of 
Tables 2-2 and 2-1 for the x >_ 2 and a _< 0 . 3 suggests that the 
distribution of T^ is converging faster to an exponential for 
Model B than for Model A. It follows from (2.2) that P{T < S n } 
is smaller for Model B than for Model A. Thus, Theorem 1 
suggests that the convergence of the distribution of T to 

X 

exponential should be faster for Model B than for Model A. 

1 -< (x) t 

4.2 The Exponential Approximation : ^ ^ e 



In this subsection simulation will be used to study the 
exponential approximation suggested by the asymptotic result 
(3.6) . This is an approximation for P{T > t} for fixed finite 

X 

x; it should be more accurate for t large. 

Two cases of Model 1 were simulated. In both cases, shock 
loads arrive according to a Poisson process with rate 1 and 
constant loads change magnitude at the times of arrival of a 
Poisson process with rate 1; the shock level magnitudes have 
an exponential distribution with mean 1. In Case A, the dis- 
tribution of the constant load magnitude is exponential with 
mean 1; in Case B, it is exponential with mean 1/2. 

In both cases considered, it is possible to determine 
analytical expressions for the integrals determining k(x) and 
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Y ( x) . The value for k(x) was found numerically. Values of 
<(x) and y(x) for the two models can be found in Tables (3-1) 
and (4-1). Note that as x increases k(x) decreases and 
approaches E[T^] . As x increases y(x) decreases and approaches 

1. As expected k(x) <_ E[T ] _1 for all levels of x. 

To assess the accuracy of the exponential approximation, 
quantiles of the simulated data were computed . These quantiles 
appear in Tables (3-2) and (4-2). For each level x, the first 
row gives the simulated quantile, the second row gives the 
approximating exponential quantile 

Q* = " F(l0" ln(Y(x) (1-a) > ; (4.3) 

and the third row gives the approximating exponential quantile 

Q* = -E[T x ] ln(l-a) . (4.4) 

The exponential approximation (4.3) is in general closer 
to the simulated quantile than (4.4). However the two approxi- 
mations become closer as x gets larger. As expected (4.3) 
approximates well the simulated quantile for a > 0.75 for 
all values of x. However if x is sufficiently large (4.3) 
approximates fairly well the simulated for a as small as 0.1. 

A comparison of Tables (3-2) and (4-2) suggests once again that 
the convergence of the distribution of T^ to exponential is 
faster for Model B than for Model A. 
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TABLE 3-1 





Values of k ( x) 


and y(x) 


for 


Model A 


Level x 


Mx) 


E [T ]~ 

X 


1 


Y (x) 


0.2 


1.708 


9 . 492 




5.942 


0.4 


1.459 


4.486 




3.254 


0.6 


1.247 


2.814 




2.364 


0.8 


1.066 


1.977 




1.923 


1.0 


0.912 


1.475 




1.670 


1.2 


0 . 780 


1.142 




1.502 


1.4 


0.667 


0.906 




1.386 


1 . 6 


0 . 570 


0.730 




1.303 


1 . 8 


0.487 


0.596 




1.238 


2.0 


0.417 


0.491 




1.193 


2.2 


0.355 


0.407 




1.155 


2.4 


0 . 303 


0 .339 




1.125 


2.6 


0.259 


0.284 




1.102 


2.8 


0.221 


0.238 




1.085 


3.0 


0.188 


0.201 




1.068 
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TABLE 4-1 



Values of k(x) and y(x) for Model B 



Level x 


k (x) 


E[T x rl 


y (x) 


0 . 2 


1.563 


4.760 


3.109 


0.4 


1.237 


2.265 


1.860 


0.6 


0.975 


1.433 


1.400 


0 . 8 


0 . 801 


1.015 


1.286 


1.0 


0.650 


0.762 


1.183 


1.2 


0.531 


0.592 


1.122 


1.4 


0 . 435 


0.470 


1.084 


1.6 


0 . 358 


0.378 


1. 059 


1.8 


0.296 


0 . 306 


1.044 


2.0 


0.244 


0.250 


1.032 


2.2 


0.199 


0.204 


1.019 


2.4 


0.166 


0.168 


1.017 


2.6 


0.136 


0.138 


1.011 


2.8 


0.112 


0.113 


1.007 


3.0 


0.092 


0.093 


1.005 
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Simulated Quantiles for Model 
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Quantile for exponential approximation ( 4 . 3 ) 
Quantile for exponential approximation (4.4) 
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